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Abstract 

We calculate the equation of state of strongly coupled Hamilto- 
nian lattice QCD at finite density by constructing a solution to the 
equation of motion corresponding to an effective Hamiltonian using 
Wilson fermions. We find that up to and beyond the chiral symmetry 
restoration density the pressure of the quark Fermi sea can be nega- 
tive indicating its mechanical instability. This result is in qualitative 
agreement with continuum models and should be verifiable by future 
numerical simulations. 
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Simulating Quantum Chromodynamics (QCD) at finite density is one 
of the outstanding problems in lattice gauge theory |]]. In fact, because 
of the sign probelm there are currently no reliable numerical simulations of 
finite density QCD with three colors even in the strong coupling limit 0. 
This is a rather frustrating situation in view of the current intense interest 
in finite density QCD fueled by the phenomenology of heavy ion collisions, 
neutron stars, early universe and color superconductivity. Therefore even a 
qualitative description of finite density lattice QCD is welcome. 

One method of studying finite density lattice QCD is to invoke the strong 
coupling approximation where analytical methods are applicable. Strongly 
coupled lattice QCD at finite quark chemical potential /i and temperature T 
has previously been studied analytically both in the Euclidean || |], ^ and in 
the Hamiltonian |6|, [7|, |8|, || formulations. One of the main objectives of these 
studies was to investigate the nature of chiral symmetry restoration at finite 
T and/or fi. This has been accomplished by constructing some effective ac- 
tion or Hamiltonian for strongly coupled lattice QCD using Kogut-Susskind 
fermions. Except for H] these effective descriptions involve composite meson 
and baryon fields which are treated in the mean field approximation.^] The 
consensus is that at zero or low T, strongly coupled lattice QCD at finite 
H undergoes a first order chiral phase transition from the broken symmetry 
phase below a critical chemical potential He to a chirally symmetric phase 
above /ic- 

In this letter we present a calculation of the equation of state of strongly 
coupled lattice QCD at finite density in the Hamiltonian formulation us- 
ing Wilson fermions. We find that up to and beyond the chiral symmetry 
restoration density the pressure of the many body system can be negative 
indicating its mechanical instability. This new result is in qualitative agree- 
ment with those obtained using continuum effective QCD models at finite 
density [[K], [T| and should be verifiable by future numerical simulations. 

As in previous studies on this subject we begin with an effective descrip- 
tion of strongly coupled lattice QCD. We shall use Smit's effective Hamil- 
tonian |T2| which involves only the quark field \I/ with a nearest neighbour 



interaction. This effective Hamiltonian has been studied in free space by Smit 
T~2f1 and by Le Yaouanc et al. []13| who subsequently extended their analy- 



sis to finite T and \x using Kogut-Susskind fermions ||. A similar effective 
Hamiltonian has recently been derived by Gregory et al. |§ to study strongly 
coupled lattice QCD at finite /i, again using Kogut-Susskind fermions. 



Henthforth we shall adopt the notation of Smit [12fl , set the lattice spacing 



1 The work of |§] does not involve composite fields but the approach is equivalent to the 
mean field approximation. 



2 



to unity and work in momentum space. Then the charge conjugation sym- 
metric form of Smit's Hamiltonian using Wilson fermions may be written 

as 



H. 



1 



off 



K 



E E« 



PlH hp4,0 



pi,...,P4 



e «((pi+P2)-n ; ) _|_ e i{(pz+PA)-ni) 



(SO;, (*t) 7 (ft)(^),(ft) - (Sl)^ (^) 5 (P3)(^) 7 (P4) 



where (Ej) = — z (707/ — ^70) with the Wilson parameter r taking on values 
between and 1. In the above Hamiltonian color, flavor and Dirac indices 
are denoted by (a, b), (a,/?) and (p, u, 7, 8), respectively, and summation 
convention is implied. N c is the number of colors. 

The three parameters in H e R are the Wilson parameter r, the current 
quark mass M and the effective coupling constant K = 2N C /(N^ — 1) 1/g 2 
where g is the QCD coupling constant. When r = M = the Hamiltonian 
possesses a U(4Nf) symmetry with Nf being the number of flavors. This 
symmetry is spontaneously broken down to U(2Nf) <g> U(2Nf) accompanied 
by the appearance of 8iV| Goldstone bosons [|12|]. A finite current quark 
mass also breaks the original U{ANj) symmetry, albeit explicitly, down to 
{2Nf) ® U(2Nf). Introduction of a finite Wilson parameter further breaks 
the latter symmetry explicitly down to U(Nf) thereby solving the fermion 
doubling problem. 

The above Hamiltonian has been derived in the temporal gauge using 
second order degenerate perturbation theory, and provides an effective de- 
scription of only the ground state of strongly coupled lattice QCD |I2|| . This 
ground state is the one in which no links are excited by the color electric 
flux. In the strong coupling limit the energy of one excited color electric flux 
link is 



E 



2N r 



1 
K 



(2) 



Therefore an extension of to finite T and/or \l will be valid as long as 
T, /j, < 1/K [Q.Q We shall see that this condition is satisfied in the present 
work. 



2 Note that in H E has been approximated by E k N c g 2 . 
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Our method for obtaining the equation of state of strongly coupled lattice 
QCD at finite density using H e ^ does not involve composite fields. Instead 
we explicitly construct a solution to the equation of motion corresponding 
to H e $ for all densities and use it to calculate the equation of state. For 



free space such a solution has been found in |L4|. This solution has the 
same structure as the free lattice Dirac field and exactly diagonalizes H c g to 
second order in field operators. It obeys the free lattice Dirac equation with 
a dynamical quark mass which is determined by solving a gap equation. 
Temporarily dropping color and flavor indices this solution is given by 

*„(t,p) = b{p)Upy~ lul{f)t + d\-p) Vu (-p)e +l ^ t (3) 

with v denoting the Dirac index. The annihilation operators for particles b 
and anti-particles d annihilate an interacting vacuum state and obey the free 
fermion anti-commutation relations. The properties of the spinors £ and rj 



are given in fl4j| . The equation of motion for a free lattice Dirac field fixes 



the excitation energy u(p) to be 

u{p) = (E sin2 (P • ni) + M 2 {p)^ (4) 

where M(p) is the dynamical quark mass. 

The extension of the method developed in [14| to finite T and \x is accom- 
plished in two steps. The first one is to make the following trivial replacement 
of the current quark mass term in H e s Eq. ([D 

Mo (To),, - M (70)^ - fioS,* (5) 

where fio is the quark chemical potential. Note that /io should not be iden- 
tified with the total chemical potential /i tot of the interacting many body 
system. As we shall see below the interaction will induce a correction to /io 
which in general is momentum dependent. We shall therefore refer to Hq as 
the "bare" quark chemical potential and treat it as a parameter. 

The second step is to observe that the annihilation operators b and d in 
Eq. @ no longer annihilate the interacting vacuum state at finite T and \i de- 
noted as \G{T,/jl)). In order to construct operators that annihilate | Q(T, //)) 
we apply a generalized thermal Bogoliubov transformation to the b and d 
operators following the formalism of thermal field dynamics |T5| 



b(p) = a p B(p) - P P B\-p) (6) 
d(p) = lp D(p) - 5 p D\-p) (7) 
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The thermal field operators B and annihilate a quasi-particle and create a 
quasi-hole at finite T and /i, respectively, while D and are the annihilation 
operator for a quasi-anti-particle and creation operator for a quasi-anti-hole, 
respectively. 

These thermal annihilation operators annihilate the interacting thermal 
vacuum state for each T and [i. 

B(p)\ Q(T, //)) = B(p)\g(T,fi)) = D(p)\g(T,fi)) = D(p)\ £(7») = (8) 

The thermal doubling of the Hilbert space accompanying the thermal Bo- 
goliubov transformation is implicit in Eq. where the vacuum state which 
is annihilated by thermal operators B, B, D and D is defined. Since we shall 
be working only in the space of quantum field operators it is not necessary 
to specify the structure of \Q(T,/j,)). 

The thermal operators also satisfy the fermion anti-commutation rela- 
tions 



(9) 



with vanishing anti-commutators for the remaining combinations. The co- 
efficients of the transformation are a n = Jl — nr, 0„ = ,/C 7„ = wl — nj 







,B(q) 


+ 


B\p) 


,B(q) 




D\p) 




+ 







and <L 



ip y 1 71^ , (3p y Tip , ^)p ^ . ,. p 

z {u p ±n)/{k B T) _|_ are p erm j distribution 

functions for particles and anti-particles. They are chosen so that the total 
particle number densities are given by 

(Q(T, fi) \ tf(p)b(p) \ Q(T, fi)) 

(g(T,n)\dft@)d(f)\g(T,n)) 



inj, where rip 



n, 



n 



(10) 
fill 



Hence in this approach temperature and chemical potential are introduced 
simultaneously through the coefficients of the thermal Bogoliubov transfor- 
mation and are treated on an equal footing. We stress that the chemical 
potential appearing in the Fermi distribution functions is the total chemical 
potential of the interacting many body system. 

In addition to these changes, we demand that our ansatz satisfies the 
equation of motion corresponding to the free lattice Dirac Hamiltonian with 
a chemical potential term given by 



p 



- ^sin(p • fii)(^i) pu + M(p)(7o) 

i 



*l(t,p),v v (t,p) 



'12) 
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As in |TJ| the mass M(p) is identified with the dynamical quark mass. Thus 
our ansatz at finite T and \x is 



a 



,B(p)-P P &(-f 



-i[cj(p)—fHot]t 



+ 



lp D\-p) - 5 p D(p) 



r] v {-p)e +i ^ + ^ t (13) 



The spinors £ and r\ obey the same properties as in free space and the exci- 
tation energy u(p) has the same form as in Eq. (|j). The unknown quantities 
in Eq. (|13|) are the dynamical quark mass and the total chemical potential. 

In this work we shall take the T — > limit which amounts to setting 
7 P = 1 and 5 P = in Eq. (^) thereby suppressing the excitation of anti-holes. 
In this limit (5 p becomes the Heaviside function f3 p = 0(fi tot — u(p)) defining 
the Fermi momentum pp where 

//tot = fe sin2 (^ ■ "0 + m2 (^») ( 14 ) 

One of the simplest quantities to calculate using the ansatz of Eq. (|13|) in 
the T — > limit is the quark number density n given by 



n 



2VN f N, 



<* 70 *> = E^tot-^(p)) (15) 



Therefore, above a sufficiently large value of /itot the quark number density 
becomes a constant which with the present normalization will equal unity. 
This saturation effect is purely a lattice artifact originating from the sm 2 (p-fii) 
term in u>(p). 

Another quantity that may be readily calculated using the T — > ansatz 
is the chiral condensate. It is found to be proprotional to the dynamical 
quark mass 

1 ■<*•> = -E«l^ us) 



2VN f N c * 1 ^ 'u(f) 

Below we shall derive a gap equation for M(p) and show that for a given 
physically reasonable set of parameters there exists a critical chemical po- 
tential above which M(p) = 0. Thus the chiral condensate may be identified 
as being the order parameter for the chiral phase transition at finite density. 

However before deriving the gap equation we shall demonstrate that in the 
T — ► limit the ansatz shown in Eq. ( |T3"D exactly diagonalizes the effective 
Hamiltonian to second order in field operators for all densities. We make use 
of the fact that our ansatz satisfies the equation of motion corresponding to 
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the free lattice Dirac Hamiltonian H° given in Eq. ([12]). Therefore we have 
the relation 



(17) 



where the symbol : : denotes normal ordering with respect to the vacuum at 
zero temperature | Q(T = 0, //)). Evaluating both sides of Eq. ( [T7| ) we obtain 

[^sin(g • n z )(7 7 i ) p5 + M(q){j Q ) pS - n tot 5 pS \ (*«»)*(*, 9 ) = 



| M (7o) p5 -/i 5 P 5 



cos (p - «r ) • d + poU 53 *)* 



+ cos (p + q ) • hA (E,)^(Si)J* + (^U^) 7 * 



Pi9 ' 



+ cos (p + q ) ■ fn ( (E,)J.(Si)^ + (E0 P ,(E0 75 



•(*««),(*, g) (18) 



with A + (p) = £(p) <g> being the positive energy projection operator 

defined in @. 

To second order in field operators the off-diagonal Hamiltonian is given 

by 



H oS \g(0,fi)) = J2) a i^M) [ M o(7o) p5 - HoSps 



p,q I 



cos(p- q) ■ nA^d^dls + (EO^CE,)^ 



+ cos (p + «f ) ■ (E^(E^ + (EO^Pi)* 



N c 4 _ _ , 

PA 1 



+ cos (p + q)-^ f (EO^CEO^ + (Ei)^(Si 
®5t (g)Dt (_^)|f?(o )A i)> 



>vs{-q) 

(19) 



From Eq. ( jl9|) we see that the elementary excitations of the effective Hamil- 
tonian are color singlet (quasi) quark-anti-quark pairs coupled to zero total 
three momentum. With the use of Eq. (|18|), the equation of motion for the r\ 
spinor and the orthonormality condition Q{p)v^{~p) 
that 



1 14] we can show 



H oS \G(0,n)) = Ma g ^(g)[-^sin(g-^)(7o7/)^ 



-M(q)(j ) vS + UtotSuS 



vs(-q 



>BU?)DU-?)\g(o,ri) 



ai q £l(q) \u{q)+fhot 



Vu(-q 



>BU?)DU-$)\g{o,ri) 







(20) 



Therefore our ansatz exactly diagonalizes the effective Hamiltonian to second 
order in field operators for all densities. 

We now derive the equations for the dynamical quark mass and the total 
chemical potential and solve them to determine our solution Eq. (|i~3|) for 
each density. To accomplish this we explicitly evaluate the right hand side of 
Eq. (|i~8"D to reveal its Dirac structure. The result may be cast in the following 
compact form 

\^sin(q ■ ni) (j y i) uS + M(q)(j ) uS - /i tot 5 uS (fy aa ) s (t,q) = 
i 

A(g)(7o 7 *U + B(q)( l0 U + C(q)5 vS ] q) (21) 



The equations for M(p) and /i to t are obtained by equating the coeffi cents of 
the 7o operator and the Kronecker delta function, respectively. 
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The gap equation determining M(p) is given by the coefficient B(q) 



M{q) = B{q 



M + ^(W)£(1-^ 



2 v y ir v v > <»(P) 



p,i 

~(1 + r 2 )cos(p + q) .n,| (22) 

The structure of this gap equation is very similar to the one in free space 
[fl 2 = 0) found in [14|. The dynamical quark mass is a constant to lowest 



order in N c but becomes momentum dependent once 1/N C correction is taken 
into account. 

Similarly, the total chemical potential is given by the coefficient C(q) 
that = ~C(q) 

= /io + ^EE^[2iV c (l + r 2 )-2(l-r 2 )cos(p + g)] (23) 
c i p 

Thus //tot is a sum of the bare chemical potential /to and an interaction 
induced chemical potential which is proportional to the effective coupling 
constant K. Furthermore, the latter contribution to //tot is momentum de- 
pendent and this dependence is a 1/N C correction just as in the case of the 
gap equation. It should be noted that the above shifting of the bare chemical 
potential by the interaction is not a new effect. For example, in the well- 
known and well-studied Nambu-Jona-Lasinio model [|16|] at finite T and /t 
the interaction induces a contribution to the total chemical potential which 
is proportional to the number density [|17], . 

The two equations ( |2"2"D and (f23|) are coupled and therefore solutions for M 
and //tot must be found self-consistently for each value of the input parameter 
/to- In Figure 1 we present M as a function of // to t for two values of K 
determined by solving Eqs. (p2|) and (|23|) to 0(N®), which is the same order 



in the 1/N C expansion used to obtain results in all previous studies of strongly 
coupled lattice QCD. At this order in N c both the dynamical mass and the 
total chemical potential are momentum independent. The values of input 
parameters are M = 0, r = 0.25 and iV c =3. 

From the figure we see that the chiral phase transition can be either first 
or second order depending on the value of the effective coupling constant. 
When K = 0.9 we find a second order phase transition with a critical chemical 



9 



CO 
CO 
CO 



CC 
Z3 

o 

o 
'E 

03 



K= 1.0 



K = 0.9 



0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 
Total Chemical Potential 



Figure 1: Dynamical quark mass M as a function of total chemical potential 
/i to t for two values of effective coupling constant K. These results were 
obtained by solving Eqs. fl22|) and (|23|) self-consistently to lowest order in 
N c using Mo = 0, r = 0.25 and N c = 3. There is a second order chiral phase 
transition when the effective coupling constant K is 0.9. The order of the 
phase transition changes to first order when K is increased to 1.0. 



potential of (/i to t)c ~ 0.716, while if the coupling constant is increased to K 
= 1.0 the phase transition becomes first order with a larger critical chemical 
potential of (/itot)c ~ 0.871. Furthermore, we find that when K = 0.9 
lattice saturation sets in around /z to t ~ 0.898 while this effect takes place 
immediately above (/i to t)c f° r K = 1.0. These values of chemical potentials 
are smaller than the energy E = 1/K required to excite one color electric 
flux link as given in Eq. (fj). Therefore with a reasonable set of parameters 
it is possible to extend Smit's effective Hamiltonian to finite density as was 
first pointed out in ||. 

Having solved for the dynamical quark mass and the total chemical po- 
tential we have constructed a solution to the equation of motion for if e ff i n 
the T — > limit to lowest order in N c . In addition we have shown that 
this solution exactly diagonalizes the effective Hamiltonian to second order 
in field operators for all densities. Therefore we may use it to evaluate the 
vacuum energy density which to lowest order in N c is given by 



i(£(0,/i)|i?eff|S(0,/i)) 



V 

-2iV c ]T< 



at 



-K(l + r 2 ) +u(p) + —rMo - //tot 
! v(p) 



+ (1 + /L 2 )/i (24) 
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Numerically we find that the difference of the vacuum energy densities in 
the symmetry restored (M = 0) and broken (M ^ 0) phases of the theory is 
positive 

^(G(P,ti) \H eS \ g(0,fi))\ M =o - ^(G(0,H) \H cS \ S(0,//))| W o > (25) 

Therefore the phase with broken chiral symmetry is the energetically favored 
phase. 

The equation of state is obtained by numerically evaluating the thermo- 
dynamic potential density using Eq. (p4|). In Figure 2 we plot the pressure as 
a function of /i tot for K = 0.9 and 1.0 with Mq = 0, r = 0.25 and iV c = 3. In 
both cases we find that the pressure of the quark Fermi sea is negative and 
monotonically decreasing in the broken symmetry phase. For K = 0.9 the 
pressure remains negative but increasing in the symmetry restored phase, at 
least until the lattice saturation point, and has a cusp where the two phases 
meet. Unfortunately, we can not make a definite quantitative statement on 
the behaviour of the pressure in the symmetry restored phase for K — 1.0 
due to lattice saturation, except to mention that there is a discontinuity 
when going from one phase to another. However, we may conclude that up 
to and beyond the chiral symmetry restoration point the quark Fermi sea 
can have negative pressure and therefore can be mechanically unstable with 
an imaginary speed of sound. 

Our conclusion regarding the (strongly coupled) quark matter stability 
at finite density is consistent with similar studies using the Nambu-Jona- 
Lasinio model ]10[ and the effective instanton induced 't Hooft interaction 



model [11]. These mean field calculations show that cold and dense quark 
matter may be unstable in the phase with spontaneously broken chiral sym- 
metry, but can become stable in the symmetry restored phase at high enough 
density. In particular, the result for the pressure of cold and dense quark 
matter obtained in JTTJ is qualitatively the same as the one shown in Fig- 



ure The possibility of unstable quark mattter lead the authors of |10| and 



TT| to speculate the formation of nucleon droplets, reminiscent of the MIT 
bag model, in the broken symmetry phase. We shall not dwell on such a 
speculation here since we are working in an artificial strong coupling regime. 
Nevertheless, it would be interesting to compare our result concerning the 
negative pressure with future lattice simulations of finite density QCD at 
strong coupling. 

In this work we studied the equation of state of strongly coupled lattice 
QCD in the Hamiltonian formulation using Wilson fermions. This was ac- 
complished by constructing a solution of the equation of motion correponding 

3 Compare Figure 1 of Q with Figure |]oi this letter. 
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Figure 2: Pressure P as a function total chemical potential // tot for two values 
of effective coupling constant K. The results were obtained using the same 
parameters as in Figure [I]. For K = 0.9, the critical chemical potential is 
(^tot)c ~ 0.716 and lattice saturation sets in around /i to t ~ 0.898, while this 
effect takes place right above (/i to t)c ~ 0.871 for K = 1.0. 

to an effective Hamiltonian which exactly diagonalizes the Hamiltonian to 
second order in field operators for all densities. We found that: the dynami- 
cal quark mass is in general momentum dependent; the interaction induces a 
momentum dependent contribution to the total chemical potential making it 
necessary to solve for the dynamical quark mass self-consistently with /i tot ; 
the elementary excitations of the theory consist of color singlet quark-anti- 
quark pairs coupled to zero total three momentum; and the broken symmetry 
phase is the energetically favored phase. 

To leading order in N c we find that the chiral phase transition can be 
either first or second order depending on the value of the effective coupling 
constant. In addition, the pressure of the strongly interacting many body 
system is found to be negative up to and beyond the chiral phase transition 
density. A similar behaviour for the pressure has been obtained with r = 
which corresponds to using Kogut-Susskind fermions. Therefore our result 
concerning the negative pressure seems to be robust, at least to leading order 
in N c , and should be verifiable by future numerical simulations of strongly 
coupled lattice QCD at finite density. 
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